Percolation of Color Sources and the Shear Viscosity of the QGP in Central A-A 

Collisions at RHIC and LHC Energies 
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The Color String Percolation Model (CSPM) is used to determine the shear viscosity to entropy 
ratio (rj/s) of the Quark-Gluon Plasma (QGP) produced in Au-Au collisions at -y/ijvjv = 200 GeV 
at RHIC and Pb-Pb at ^/snn = 2.76 TeV at LHC. The relativistic kinetic theory relation for 
rj/s is evaluated using CSPM values for the temperature and the mean free path of the QGP 
constituents. The experimental charged hadron transverse momentum spectrum is used to determine 
the percolation density parameter £ in Au-Au collisions (STAR). For Pb-Pb at ^/sjvjv = 2.76 TeV 
£ values are obtained from the extrapolation at RHIC energy. The value of rj/s is 0.204±0.020 
and 0.262±0.026 at the CSPM initial temperatures of 193. 6±3 MeV (RHIC) and 262.2 ±13 MeV 
(LHC) respectively. These values are 2.5 and 3.3 times the AdS/CFT conjectured lower bound 1/4-k. 
We compare the CSPM rj/s analytic expression with weak coupling (wQGP) and strong coupling 
(sQGP) calculations. This indicates that the QGP is a strongly coupled fluid in the phase transition 
region. 

PACS numbers: 12.38.Mh; 25.75.Nq 
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The observation of the large elliptic flow at RHIC in 
non-central heavy ion collisions suggest that the matter 
created is a nearly perfect fluid with very low shear vis- 
cosity [3-[l] . Recently, attention has been focused on the 
shear viscosity to entropy density ratio rj/s as a mea- 
sure of the fluidity p-Q- The observed rj/s based on 
viscous hydrodynamics analyses at RHIC, is suggestive 
of a strongly coupled plasma [H, [l(| • The effect of bulk 
viscosity is expected to be negligible. It has been conjec- 
ture, based on infinitely coupled supersymmetric Yang- 
Mills (SYM) gauge theory using the correspondence be- 
tween Anti de-Sitter(AdS) space and conformal field the- 
ory (CFT), that the lower bound for rj/s is 1/47T and is 
the universal minimal viscosity to entropy ratio even for 
QCD [11]. However, there are theories in which this lower 
bound can be violated [l2j . 

In this letter, we explore the color string percolation 
model (CSPM) [II HI] to obtain rj/s as a function of the 
temperature above and below the hadron to QGP transi- 
tion. The experimental values are for Au-Au collisions at 
y/sl^i = 200 GeV at the RHIC and for Pb-Pb collisions 
at ^/JWn = 2.76 TeV at the LHC. 

The CSPM describes the initial collision of two heavy 
ions in terms of a dense partons of interacting colored 
strings which start to overlap to form clusters [l3|, [l4| ■ 
The interactions of strings reduces the hadron multiplic- 
ity ji and increases the average transverse momentum 
\J (p 2 ) of these hadrons, so that the total transverse mo- 
mentum is conserved. The observables fi, and {pf) are 
given by 
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where (Iq is the multiplicity of hadrons from a single 
string and (pf)i is the average transverse momentum 
squared of a single string. N s is the number of longi- 
tudinal strings formed in the collisions and F(£) is the 
color suppression factor [l| [l4| 
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We note that CSPM is a saturation model where 
(jpf)i/F(^) is similar to the saturation momentum scale 
Ql in the Color Glass Condensate (CGC) model [l5j . 

In two dimensional percolation theory the relevant 
quantity is the dimensionless percolation density param- 
eter given by [HI [l3| 
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where 5*i is the transverse area of the single string and 
Sn is the transverse nuclear overlap area. The critical 
cluster which spans Sn, appears for £ c > 1.2 [16j. To 
evaluate the initial value of £ from data for Au-Au colli- 
sions, a parameterization of pp events at 200 GeV is used 
to compute the pt distribution fllHl~9| . 



dN c /dp 2 t =a/(po+Pt) c 
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where a, po, and a are parameters used to fit the data. 
This parameterization also can be used for nucleus- 
nucleus collisions to take into account the interactions 
of the strings [l4[ 
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where S n corresponds to the area occupied by the n over- 
lapping strings. The thermodynamic limit, i.e. letting n 
and S n — > oo while keeping £ fixed, is used to evaluate 



e~*)SN- Thus the effective no. of sources is given by the 
total area occupied by the strings divided by the area of 
an effective string SiF(£). 



dN c /dp 2 t 
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In pp collisions (nSi/ S n ) pp ~ 1 at not very high energies 
due to the low overlap probability. 

The connection between £ and the temperature T(£) 
involves the Schwinger mechanism (SM) for particle pro- 
duction. The Schwinger distribution for massless parti- 
cles is expressed in terms of p\ 
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with the average value of string tension, (x 2 ). Gaussian 
fluctuations in the string tension around its mean value 
transforms SM into the thermal distribution 120] 



dn I dp 2 

with (x 2 ) =7T(p 2 ) 1 /F(0. 

The temperature is expressed as []J 
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The string percolation density parameter £ which charac- 
terizes the percolation clusters also determines the tem- 
perature of the system. In this way at £ c = 1.2 the perco- 
lation phase transition at T(£ c ) models the thermal de- 
confinement transition. We will adopt the point of view 
that the experimentally determined chemical freeze-out 
temperature is a good measure of the phase transition 
temperature, T c (22|. (p 2 )i is calculated at £ c = 1.2 using 
the universal chemical freeze-out temperature of 167.7 ± 
2.6 MeV H|. This gives y/(p 2 )i = 207.2 ± 3.3 MeV 
which is close to ~ 200 MeV used in a previous calcula- 
tion of the percolation transition temperature plj . 

The relativistic kinetic theory relation for the shear 
viscosity over entropy density ratio, r]/s is given by [9Tl24j 
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where T is the temperature and A m / P is the mean free 
path, n is the number density of an ideal gas of quarks 
and gluons and Otr the transport cross section. 

In CSPM the number density is give n by the effective 
number of sources per unit volume [291 ] 
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L is the the longitudinal extension of the string L ~ 1 
fm. The area occupied by the strings is related to the 
percolation density parameter £ through the relation (1 — 
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In general N sources is smaller than the number of single 
strings. N sources equals the number of strings N s in the 
limit of £ = 0. The number density of sources from Eqs. 
(12) and (13) becomes 



(14) 



In CSPM the transport cross section a tr is the trans- 
verse area of the effective string S\F(£)- Thus a tr is di- 
rectly proportional to , which is in agreement with the 
estimated dependence of a tr on the temperature [24], [2|| . 
Equations (11)-(15) give 
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FIG. 1: Color suppression factor i ? (£) as a function of 
^j^/SN(fm~ 2 ). The solid red circles are for Au-Au colli- 
sions at 200 GeV(STAR data) [TJ]. The error is smaller than 
the size of the symbol. The line is fit to the STAR data. The 
solid blue squares are for Pb-Pb at 2.76 TeV. 

In the above expression for rj/s the key quantity £ is 
obtained from the determination of F(£) given in Eq. 
(7). 

Figure 1 shows a plot of as a function of charged 
particle multiplicity per unit transverse area /Sn 
for Au+Au collisions at 200 GeV for various centralities 
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FIG. 2: Energy density e as a function of percolation density 
parameter £. The line is fit to Au-Au data. The LHC point 
is obtained from e/T 4 value [T^ |. 



TABLE I: The percolation density parameter £, temperature 
T, the energy density e and the ratio of shear viscosity to 
entropy rj/s for percolation threshold ( £c ), Au-Au at 200 
GeV and Pb-Pb at 2.76 TeV (estimated) [l3,[3- Au-Au is 
for 0-10% and Pb-Pb is for 0-5% central events. 



System 




T (MeV) 


e(GeV/fm a ) 


rj/s 


Critical 


1.2 


167.0 ±2.6 


0.94±0.07 


0.210±0.016 


Au-Au 


2.88±0.09 


193.6 ±3.0 


2.27±0.16 


0.204±0.020 


Pb-Pb 


10.56±1.05 


262.2±13.0 


8.32±0.83 


0.260±0.026 



the Equation of State (EOS) it was demonstrated that 
£ is proportional to energy density e at RHIC energies 
17]. Figure 2 shows a plot of e vs £ for Au-Au collisions 
at Jsnn = 19.6, 62 and 200 GeV for 0-10% centrality 
[l7l Il9l. 27 1 . The energy density for the LHC energy was 
calculated using the value of e/T 4 = 12.5 at T = 262.2 
MeV [13 ■ This value is shown in Fig. 2. The line is fit to 
the STAR data. This shows that the linear relationship 
between £ and e is valid even at 2.76 TeV when extrapo- 
lated from RHIC energy. Table I gives the CSPM values 
for £, T, e and r]/s at T/T c =1, 1.16 and 1.57. 

The CSPM r\j s values are obtained using Eq.(15). This 
ratio grows slowly as £ 1 / 4 and remains small between 
RHIC and LHC energies. Figure 3 shows a plot of rj/s 
as a function of T/T c calculated for different values of 
£. The experimental value for Au-Au at 200 GeV for 
0-10% centrality is marked with arrow at T/T c = 1.15. 
The estimated value of rj/s for Pb-Pb is also shown in 
Fig. 3 at T/T c = 1.57. The lower bound is shown in 
Fig.3 given by AdS/CFT [lj. The results from STAR 
and ALICE show that the rj/s value is 2.5 and 3.3 times 
the KSS bound [Il|. The small rj/s in CSPM is due to 
the small transverse correlation length, which decreases 
as the energy increases compensating the increase of the 
effective number of sources in such a way that the X m f p 
given in Eq. (11) remains almost constant [29| . 

The viscosity to entropy ratio has been calculated for 
both the weakly (wQGP) and strongly (sQGP) coupled 
QCD plasma Q. In this letter the CSPM results are 
compared with these calculations of wQGP and sQGP 
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[HHIi. The error on F(£) is - 3%. F(£) decreases in 
going from peripheral to central collisions. The £ value 
is obtained using Eq. (3), which increases with the in- 
crease in centrality. The fit to the Au-Au points has the 
functional form 



where 



F(£) = exp[-0.165 
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Recently, ALICE experiment at LHC has published the 
charged-particle multiplicity density data as a function of 
centrality in Pb-Pb collisions at ^/snn = 2.76 TeV [26|. 
The STAR results for Au-Au collisions at ^/snn = 200 
GeV have been used to estimate F(£) values for Pb-Pb 
collisions at different ccntralities using the fit function 
Eq.(16) for Au±Au. The ALICE data points are shown 
in Fig.l. For central 0-5% in Pb-Pb collisions £ = 10.56 
as compared to £ = 2.88 for Au-Au collisions at 200 GeV. 
The temperature was obtained using Eq. (10). For Pb- 
Pb collisions the temperature is ~ 262.2 MeV for 0-5% 
centrality, which is ~ 35 % higher than the tempera- 
ture from Au-Au collisions fl7j . In our previous work on 
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In the sQGP phase the ratio is 
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Both sQGP and wQGP rj/s values are shown in Fig. 
3. It is seen that at the RHIC top energy r]/s is close 
to the sQGP. Even at LHC energy it follows the trend 
of the sQGP. Extrapolating the rj/s CSPM values it is 
clear that it will approach the weak coupling limit near 
T/T c ~ 5.8 as the temperature is increased. Figure 3 
shows that both Au-Au and Pb-Pb systems are close to 
the strongly interacting QGP calculations @. 

The CSPM model calculations have successfully de- 
scribed the dependence of the nuclear modification fac- 
tor on the azimuthal angle and obtain a similar pt de- 
pendence of v2 at RHIC and LHC energies 0, [2{| ■ In 
addition CSPM coupled to a Bjorken hydrodynamic ex- 
pansion determined the equation of state of the QGP 
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FIG. 3: r)/s as a function of T/T c . Au-Au at 200 GeV for 
0-10% centrality is shown as solid black square. wQGP and 
sQGP values are shown as dotted blue and green lines respec- 
tively [jj. The estimated value for Pb-Pb at 2.76 TeV for 
0-5% centrality is shown as soild blue square. The red dotted 
line represents the extrapolation to higher temperatures from 
the CSPM. 



and the bulk thermodynamic value of e/T 4 and s/T 3 
which are in good agreement with Lattice QCD simula- 
tions [H El ■ 

In summary the relativistic kinetic theory relation for 
shear viscosity to entropy density ratio rj/s = X m fp 
was evaluated using the Color String Percolation Model. 
The mean free path X m f p is proportional to the number 
density of the colored sources in CSPM. The color sup- 
pression factor -F(£) was extracted from the transverse 
momentum spectrum of charged hadrons. We found rj/s 
= 0.204 ±0.020 at T/T c = 1.15 ( RHIC ) and r//s =0.260 
±0.020 at T/T c = 1.57 (LHC). The CSPM tj/s was com- 
pared to weak coupling (wQGP) and strong coupling 
(sQGP) calculations. The CSPM values are close to the 
sQGP estimates at RHIC and LHC energies but could 
approach the wQGP estimate near T/T c ~ 5.8. In the 
phase transition region rj/s is 2-3 times the conjectured 
quantum limit for RHIC to LHC energies. 
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